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A general theory of optical forces on moving bodies is here developed in terms of generalized/4 x 4 
transfer and scattering matrices. Results are presented for a planar dielectric multilayer of arbitrary 
refractive index placed in an otherwise empty space and moving parallel and perpendicular to the 
slab-vacuum interface. In both regimes of motion the resulting force comprises lateral and normal 
velocity-dependent components which may depend in a subtle way on the Doppler effect and TE- 
TM polarization mixing. For lateral displacements in particular, polarization mixing, which is here 
interpreted as an effective magneto-electric effect due to the reduced symmetry induced by the 
motion of the slab, gives rise to a velocity dependent force contribution that is sensitive to the phase 
difference between the two polarizations amplitudes. This term gives rise to a rather peculiar optical 
response on the moving body and specific cases are illustrated for incident radiation of arbitrarily 
directed linear polarization. The additional force due to polarization mixing may cancel to first 
order in V/ c with the first order Doppler contribution yielding an overall vanishing of the velocity- 
dependent component of the force on the body. The above findings bare some relevance to modern 
developments of nano-optomechanics as well as to the problem of a frictional component to the 
Casimir force. 
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I. INTRODUCTION 



Although Minkowski derived the basic equations for 
describing the electrodynamics of moving media in 
1908 [l], it is clear that there is still a good deal to be 
done to understand the effect of the electromagnetic field 
on the motion of a dielectric body. While the Abraham- 
Minkowskii debate [2H3] is the most persistent instance of 
this, a more tractable problem can be found within the re- 
cent controversy surrounding the existence of a frictional 
component to the Casimir force [5HPT] . Although a quan- 
tum phenomenon, much of the initial debate was con- 
cerned with classical radiation pressure. Namely, what 
is the electromagnetic force on a dielectric medium in 
motion? The concern was with the relative importance 
of the relativistic transformation of polarization, versus 
that of frequency (the Doppler shift) to the force exerted 
by the electromagnetic field on a moving medium 8 10 . 

A further motivation for the current investigation 
comes from the remarkable recent progress in optomc- 
chanics (see [H]) which uses the motional effects of ra- 
diation pressure to cool mechanical degrees of freedom 
to sub-Kelvin temperatures |13H15j . and even close to 
the quantum ground state [IB]. Relevant to this is work 
on 'radiation damping' jTTJ, [TS] in highly dispersive sys- 
tems such as photonic crystals [TS], and atomic multi- 
layers [20 . This paper adds to other recent work on 
velocity dependent forces in optomechanics |21) . gener- 
alizing such treatments to an arbitrary direction of mo 
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tion. It is also worth mentioning that our findings may 
be applicable within the field of metamaterials (see for 
example [22 -24 ), where some aspects of motion can be 
mimicked with a suitably engineered magnetoelectric re- 
sponse (however, see |25j). and where optical proper- 
ties such as extreme sensitivity to polarization (see sec- 
tion 



III A ) may be achieved with structured surfaces 



Here we focus on the theory of radiation pressure for 
planar media in relative motion, and show how the mo- 
tion leads to velocity dependent forces that always have 
two distinct origins; one contribution coming from the 
mixing of s and p polarizations by the motion; and 
the other more familiar effect coming from the Doppler 
shifted frequency that appears within the rest frame 
reflection and transmission coefficients. We find that 
whereas the latter effect may be amplified through an in- 
crease in the dispersion of the medium, the former may 
be amplified through an increase in the sensitivity of the 
medium to the polarization at small angles of incidence. 
We quantify the relative importance of these two effects 
and show situations where either contribution can be sig- 
nificant. 

In section [TT] we outline the general theory of radiation 
pressure for planar media, where the response is charac- 
terized in terms of scattering, or transfer matrices. Sec- 
tion |III| then applies this to the case of moving media, 
illustrating that the motion can be understood in terms 
of a linear transformation of the scattering matrix. We 
comment briefly on the likely implications of our findings 
to the vacuum force experienced by two dielectric plates 
in lateral motion. 
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II. CALCULATING RADIATION PRESSURE 
USING THE TRANSFER MATRIX 

Optical forces on a body are assessed by calculating 
either the net Lorentz force experienced by the body, 

~ = J[pE+jxB]d 3 x, (1) 

where p = -V-Pand j = VxM+dP/dt (c.f. H3HH]), 
or the rate of momentum being lost from the field in the 
space outside the body. The rate of momentum leaving 
the electromagnetic field can be written as an integral 
over the surface of the energy-momentum tensor (29H31j , 

rpflV _ ( ®F \ 

~ \cV F -cr F J 

= ff (E 2 +c 2 B 2 ) ce E xB \ 

\ ce ExB t^ F -e Q [E®E + c 2 B®B\) 

(2) 



where S v and "P F are the energy and momentum den- 
sity, respectively, and erp is the Maxwell stress tensor. 
In this paper we choose to work in terms of the energy- 
momentum tensor, for it allows the optical force to be 
straightforwardly written in terms of the medium reflec- 
tion and transmission coefficients. Note that throughout 
it is assumed that the region outside of the material has 
e ~ eo, and p ~ /io, so that we may avoid possible ambi- 
guities in the form of the energy momentum tensor within 
the surrounding space [2j 0] . 



A. The energy-momentum tensor for a plane wave 
reflecting from a surface. 

For a field with a harmonic time dependence 
E ul (x)e~ zu ' t and £? w (x)e _J "*, observed over a time pe- 
riod At, where 1/uj -C At, we can deal with the time 
average of ^ over many optical cycles. Using the fact 
that only the real values of the fields enter the energy- 
momentum tensor, the time average of p) IS, 



Re 



£11 

4 



| 2 + c 2 |£ 
\E U x B* 



2 



[E u x Bt 



(T 00 )l 3 - f [E u ® El + c'B^ ® Bt] 



(3) 



For the specific case of a plane wave reflecting off a sur- 
face, as shown e.g. on the left hand side of the planar 
slab in fig. I, the electric and magnetic fields are, 

^w = -^EE(- 1 ) 9 4 ±) 4 ±) ^ <±, ' K (4) 

where the wave-vectors 



fc (±) = ±\k x \x + k\\k\\ 



(5) 



are decomposed as usual in terms of their normal and 
in-plane components. For each wavevector there are two 
transverse polarization directions labelled by q = 1,2 
with q — (q mod 2) + 1. These polarization vectors are 
defined as, 

= x X fey = [k y z — k z y] s-polarization (6) 
normal to the plane of incidence (TE polarization), and 



k 



C±) 



X e\ = ClI^T Cxfcll p-polarization (7) 



Here ot\ 2 are the electric field complex amplitudes of 
the right (+) and left (— ) propagating waves, of either 
s (TE), or p (TM) polarization. The summation in Q 
is over the modes propagating to the right (+) and to 
the left (— ), and over the two polarizations labelled by q. 
The geometrical details of the slab and the plane-wave 
configurations at the two interfaces on the left and the 
right are shown in figure [TJ where the origin of the x axis 
is taken at the boundary between the medium and the 
free-space on the left-hand side of the slab and the y and 
z axis lie parallel to the interface. 

We proceed to use ^ to calculate the radiation pres- 
sure experienced by the slab from a plane wave reflecting 
on the left. Using the form of the fields given in Q, 
we can write down the general energy-momentum tensor 
outside a planar surface in terms of the field amplitudes, 
a 1*2 ■ T ne time-averaged energy density from |s| is found 
to be, 



£0 
2 



E 



El 
{±> 



(9) 



parallel to the plane of incidence (TM polarization) , with 
the directional terms 



Cx — c\k x \/uj = cos 9 Gi = ckn/uj = sin( 



(8) 



where the terms A qr {x) = 2^{a ( q +) a[ > e 2 ^ x ) are in- 
terpreted as being due to a stationary interference pat- 
tern between incoming and outgoing waves. The time- 
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FIG. 1: Geometry of the slab interfaces showing the input 
(red) field amplitudes a4 and ot^' and the output (blue) 
amplitudes £>£ and \ The slab surface normal is along 
±x while the orientation of the two polarizations (1, 2) are in- 
dicated in the input amplitudes. The inset shows the unit po- 
larization vectors |6j|7f and their relative orientation in terms 
of the angle defined in (ISj) . 



averaged momentum density is similarly, 



<w = |E 

{9> 



E 

{±> 



v(±) 



2 fe (±) 



C|| A,,(a;)fe|| +(-l)' 1 C x A qq (x)e 1 



(10) 



where it is clear that there is a momentum flux parallel 
to the surface as well as normal to it. The calculation 
of the 3x3 tensor, —(ap), is slightly more involved, yet 
noting that 



e\ (x> e\ + e 



— 13 = — fc (g) fc , 



The time-averaged energy-momentum tensor given 
above allows us to determine the force exerted by the 
incident plane-wave on the medium, as well as the 
work done on it. We consider fields in the rest frame 
of the medium and calculate the time-averaged four- 
momentum lost from the field per unit time in the space 
outside of the material. The rate of change of the four- 
momentum density V Y = (<S Y /c. "P F ) at any point outside 
of the medium is related to the spatial divergence of the 
energy-momentum tensor, T piy , via the vanishing four- 
divergence, dfjT^ 11 = 0, i.e, 



1 dT 0u dVl 



dt 



dt 



dx l 



(12) 



Integrating ( 12 ) over the volume of space outside of the 
material, V, gives us the rate of change of the total four- 
momentum, P v , in free space, 



dt 



Qrpxv 
dx % 



d A x = - 



T iv dS u 



av 



where the three dimensional divergence theorem has 
been applied to obtain the final expression. Due to the 
global conservation of four-momentum, whatever is lost 
in the free space is taken up by the region of space occu- 
pied by the medium, 



dt 



dPl 
dt 



T w dS l 



al- 



and, 



Upon averaging over a time interval of many optical cy- 
cles this change of momentum is that taken up by the 
dielectric medium (the rate of change of electromagnetic 
momentum within this region of space averages to zero 
over such a time period). Writing out the components 
explicitly using the notation introduced in ([2]), we have, 



ex <g> ei 



one has. 
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(11) 



Equations ([9 11) determine the energy-momentum ten- 
sor for the case of a plane wave reflecting from the planar 
surface of a generic medium at rest, taking into account 
both polarizations. 



dE u 
~dt 
dP M 
dt 



c- 1 (V F ) ■ dS 
'av 



■dS, 



av 



(13) 



where each surface element, dS points into the material 
(c.f. figure [lj. In general, for an object of finite size in 
a plane wave field, there will be reflection from each sur- 
face of the object that will complicate the evaluation of 
the reflected fields and of the integrals in (13 1. We as- 



sume throughout the case of a thin planar slab sufficiently 
extended over the y — z plane so that translational sym- 
metry along this plane is preserved and the total force 
can be calculated integrating only over the left and right 
interfaces parallel to this plane. 



With these assumptions, and an application of ( 10 ) 



and (11) to (13), the components of the time-averaged 



four-force experienced by such a medium in its rest frame 
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are, 



and the transfer matrix has been compactly written as, 



dE M 
dt 

dPu 

dt 



Ae ( x 



EE± 

{«> {±> 

EE± 

{<?} {±} 



y {±) 



'(±) 



(14) 



where A is the integration area (the cross section of the 
beam) . The subscripts L and R indicate here the electric 
field amplitudes (a) on the left (x < 0), and right (x > 
d) of the slab, respectively (figure [l]). Hence, the time 
averaged four force on the slab is determined once the 
field amplitudes are known in vacuum. For a given plane 
wave of frequency ui and wavevector fc,the two results in 



( 14 ) may be combined to yield the rather sound relation, 



dP* 
dt 



k I dE 



dt 



which holds true for either polarization. 



B. 4x4 Transfer and Scattering matrices. 



a 



In gen eral , the relations between the amplitudes, 
tfL/fl * n jl^j ma y b e obtained by employing transfer or 
scattering matrix methods (e.g. |31H33| ). which we now 
use — along with the input fields — to determine the radi- 
ation pressure experienced by a general dielectric slab. 
The transfer matrix approach is ideally suited for the 
consideration of media inhomogeneous in one direction 
only (here taken along x) as the medium can be split 
up into as many layers as accuracy demands, and the 
transfer matrices of each of these layers applied one after 
the other to obtain the matrix associated with the whole 
slab. We will be using a 4 x 4 transfer matrix theory as is 
commonly applied to anisotropic media (e.g. [31]). This 
is because even if a medium is isotropic in its rest frame, 
when set in motion its symmetry is lowered. 



1. The transfer matrix 

A typical 4x4 transfer matrix, T, connects the field on 
the left of a slab (x < 0) to the field on the right (x > d 
see figure [IJ and is such that, 



OtlR 
OC2R 



Tn Ti 2 \ ( ol\l 

T 2 l T22) \Ot2L 



(16) 



where the field amplitudes have been written in terms of, 



*qL/R 




Tn T12 
T21 T22 



ftll ti2 tl3 tu\ 

t21 t22 ^23 ^24 

t31 £32 ^33 ^34 

V41 £42 t 43 £44/ 



(17) 



with each sub-matrix, T qr relating amplitudes with po- 
larization r on the left hand side of the medium, to am- 
plitudes with polarization q on the right hand side. 
The elements of T in (17) are not all independent 



but in general depend on specific symmetries that the 
medium may posses. Most importantly, symmetries that 
hold in the medium rest frame may be lifted when the 
medium is set in motion — three specific cases of this are 
shown below. 

When the plane of incidence is a mirror symmetry 
plane, then s (TE) and p (TM) polarizations do not mix, 
i.e. are decoupled and the off diagonal sub-matrices of 
(fl7|) will vanish, T i2 = T 2 i — 0, reducing the transfer 



(15) matrix in (16) to the direct sum: 



Tu 



22 



Tu 
T 22 



This is e.g., the usual case for an isotropic medium at 
rest, for which a 2 x 2 transfer matrix can instead be 
used, but no longer applies in general when the medium 
is set into motion. The plane of incidence remains in fact 
a plane of mirror symmetry for a medium moving along 
x, but not for a medium moving along y (unless k z = 0) 
in which case the polarizations mix with non vanishing 
off diagonal elements in (17), as will be discussed below. 



Meanwhile for a magneto-electric medium at rest, the 
constitutive relations take the form D = e • E + r\ • B 
and H = fjT x • B + £ • E |30l EH [36] and, in general, 



all entries in (17) are non-zero. We may therefore refer 



to the polarization mixing induced by motion as to a 
"magnetoelectric effect" [51 130j. 

The time reversal for a plane wave is obtained by 
changing fen into —fey and into (a q T ^)*. When 

time reversal symmetry applies, e.g. for a lossless non- 
magnetized medium at rest, the transfer matrix should 
be the same when waves are interchanged in this way, 
i.e. left-going waves viewed with time running in re- 
verse should be reflected and transmitted as if they were 
right-going incident waves with time running forwards, 
provided that changing the sign of k y and k z is also im- 
material as for an isotropic medium. Performing such a 
transformation on the transfer matrix leads to, 



& xT u(T x 

xT '21" ' X 



rrrk 

v- 1 12 



cr xT2 2 cr o 



Tu T12 
T21 T22 



(18) 



where cr x is the usual Pauli matrix (er^ = f 2 ). There- 
fore, ( 18 ) shows that time reversal symmetry requires 
that the sub-matrices of T transform as, Tjj(feii) = 
a x T*A—k\\)a x . For example, 
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in 

*21 



*12 
f 2 2 



'22 

t* 
'"12 



t* 
'21 

t* 

Hi 
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For a moving medium, time reversal symmetry does not 
apply as it would require the velocity of the medium to 
change sign. 

If the median plane of a slab is a plane of mirror sym- 
metry, an input of a given amplitude produces the same 
reflected and transmitted amplitudes whether incident 
from the right or the left of the slab, and we refer to this 
as left-right symmetry. Noting that T transfers ampli- 
tudes on the left of the slab to the right, and that T~ 
transfers amplitudes from right to left, this condition is 
equivalent to, 



T 21 



T22 











T 2 i 



T i2 
T22 



a x 
er, 



(19) 

where the Tij are the sub-matrices of the inverse, T . 
The cr x matrices are present to interchange left and right 
propagating waves, as we should in examining the mirror 
symmetric situation. The reason for the minus sign in 
one of these matrices (— <r x ) in (191 is more subtle: it is 
present due to the definition of p-polarization |7]) (see 
figure [T]) . If the polarization of the incoming p-polarized 
wave on the right hand side of the slab points away from 
the surface then that on the left hand side points into the 
surface. Therefore in examining the mirror symmetric 
situation we should multiply the p-polarized amplitudes 
by -1. 

In general, the conditions for left-right symmetry aris- 
ing from ( 19 ) are complicated as they require the calcula- 
tion of the inverse of a 4 x 4 matrix. However, in the case 
of media for which T — TnQ)T 22, the condition becomes 
rather compact. In this case it requires the elements of 
Tn to obey, 



fer matrix to an experimental situation. In general the 
known quantities are the 'input' field amplitudes, 

and a„ R , rather than either of the quantities, a q L/R 
(see figure [I]) . A scattering matrix, S, rather than trans- 
fer matrix relates input and output, 



a l (OUT) 
Q!2 (OUT) 



Su S\2 
$21 S22 



OL\ (in) 
a 2 (IN) 



(21) 



where, 



Su S12 
S21 S22 



(Tn Tin T12 K u \ 

Tin Tn 7li2 Tn 

T.1 7Z-21 T22 7^-22 

\R-21 Ti H22 T22 / 



(22) 



with TZ qr (TZ qr ), Tqr(Tq r ) representing reflection and 
transmission coefficients for g-polarized waves, given an 
r-polarized wave impinging from the left (right), and 
where, 



a 



"■q (IN) 



(23) 



and, 



l q (OUT) 



l qR 



(24) 



Tn 



tn ti2 

til ^22 



det(T 



in 

-*12 



-*21 
t22 



(20) 



with an identical relation also holding for the matrix T 2 2- 
For this particular case, the Tn matrices must be anti- 
symmetric and have a determinant equal to one (when 
such a medium is also time reversible, i 2 2 = and 
t2i = t* 2 = —ti2)- For radiation passing from vacuum 
into a such a medium and then back into vacuum (as con- 
sidered here), this determinant will always equal unity. 
We finally note that when a medium is moving along x, 
left-right symmetry is lifted, while it still applies if the 
motion is perpendicular to x. 



2. The scattering matrix 

Although undoubtedly useful for calculations, there 
is a slight awkwardness in the application of the trans- 



Direct evaluation of the scattering matrix for a generic 
slab is in general rather involved. One would then pro- 
ceed, first finding the transfer matrix for the whole slab, 
and then converting this into a scattering matrix for the 
calculation of the force, via (14). 



The relationship between the elements of the scattering 
matrix ( 22 ) and the transfer matrix can be found through 



expanding the matrix equation ( 16 ) in terms of the 16 



elements of ( 17), and writing the output field amplitudes, 
(24), in terms of the input amplitudes and the tij. This 



gives, 



r(") 



r <+) . T 



)•(• 



.(+) 



r(") 



where we have introduced two matrices, 7r (+) = 
diag(l, 0, 1, 0) and = diag(0, 1, 0, 1). Explicitly eval- 
uating the components of the S matrix we have 
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(^22^44 — ^24^42) 



I M33 (tl2*44 — il4*42) — M31 (^22*14 — ^12i24) N 

(^24^41 — ^44^21) ^44 (^24^43 — ^44^23) ~^24 

— Mi 3 (^32^44 — ^34^42) -^11 (^22*34 — ^32*24) 

\(^42*21 — ^22*4l) ~ ^42 (^42^23 — ^22^43) ^22 



(25) 



where is the minor of T, e.g. 



M 
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^21 ^22 ^24 
^31 £32 £34 
£41 £42 £44 



Note that, as expected, when the magnetoelectric effect is 
not present, the scattering matrix also reduces to a direct 
sum, S — 5n®S l 22, with the comparatively simple form, 



S n 



1 

t 2 2 



1 

"*21 



*12 
f 



'22 



1 

£44 



1 t 3 4 
"*43 1 



where we have assumed that det(Tn) = det(T 2 2) 
I . From (251 it is clear that in the general case there 



is quite an intricate relationship between the scattering 
and transfer matrices. We also note that in the case of 
systems with a resonance in reflection or transmission 
(e.g. a Fabry-Perot cavity), the resonant frequency is 
generally determined by £22^44 = ^24^42, rather than the 
usual separate conditions for the two polarizations, £22 = 
and £44 = 0. In addition, if the medium is left-right 
symmetric, one has t 2 i — —t\2 and £43 = — £34, so that; 
72n = 72-11 = inl^Tij 71 1 = 71 1 = I/E22; 72.22 = 72-22 = 
£34/^44; and, T22 = T22 = 1/^44- Due to the restrictions 
of the previous section (ti2 = — t\^), time reversibility 
has the further consequence that, TZ qq T qq + T qq lZ qq = 0. 

Finally, in terms of the scattering matrix elements de- 
fined in (21 1 the four-force (14) can be rewritten as, 



dE u 
~dt 

dP M 

dt 



AcqCx 



( IN ) 



1 14 - S^SJ ■ a (1N) 



(26) 



where the matrices in the second equation operate on the 
input amplitudes outside of the vector, and R = er z ©er z , 
with <j z the usual Pauli matrix. Hence, once we have 
determined the transfer matrix for any given planar slab 
we have, via (25-26), thereby determined the four-force, 
for any given input field. In ( 26 ) only the input field 



amplitudes enter, which is the main advantage of ( 26 ) 



over (14). When the slab is lossless, the scattering matrix 



(22), to give the usual result (e.g. within [37] these ex- 
pressions appear for a single polarization in the force ex- 
perienced by a transparent dielectric slab), 



dE M 
~dT 

dPu 
dt 
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1 - 172 



2 -\T q 12 



'991 
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[1 + \K q 

9L| Ui-irc* 



<+) 



\T qq \ 

\Taa\ 



C||fe|| 



(27) 



It is clear that the lateral force is simply proportional 

1 1 2 1 1 2 1 1 2 

to the absorption \A qq \ = (1 — \7t qq \ — \T qq \ ), while 

the normal force contains, besides the contribution of the 

absorption, a term proportional to twice the reflectivity, 

i.e. 2 \H qq \ + \A qq \ = (1 + \H qq \ - \T qq \ ), as reflec- 
tion entails the exchange of twice the momentum in the 
direction normal to the surface. 

To conclude this section we note that in the above we 
have not yet assumed that the dielectric is in motion. 
We shall find in the following section that for certain 
regimes, the motion of the medium induces a coupling 
between the s and p polarizations. However, the above 
conclusions hold for any media that couple polarizations, 
moving or otherwise. 



III. MOVING MEDIA AND RADIATION 
PRESSURE 

In the following we use the formalism of section |H] to 
study radiation pressure forces on a moving slab, and we 
will specifically examine two regimes of motion, namely, 
when the slab moves parallel (A.) and perpendicular (B.) 
to the interface. 



A. Motion in the y—z plane 

When an incident plane wave of a frequency oj reflects 
from a surface that moves in the plane orthogonal to 
the surface normal, the reflected plane wave also has fre- 
quency oj, and the field is monochromatic in both rest 
frame and lab frame. Therefore we may apply the re- 



is unitary, S = I4, and any lateral force on the slab 
will vanish. 

For the simplest case of incidence from the left (a^ H — suit ( 26 ) directly. Without loss of generality, we take the 
0), onto a reciprocal, non-magneto-electric medium, 
equation (26) can be expanded, using the notation of 



motion to be along y to be representative of the general 
motion in the y-z plane (see figure^]). 
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FIG. 2: For certain angles of incidence, a laterally moving 
medium mixes polarizations in transmission and reflection. 
For instance, an s-polarized input will generate p-polarized 
reflected and transmitted amplitudes so long as k z 7^ 0. In all 
cases the Doppler shift makes the reflection and transmission 
coefficients dependent upon k y in such a way that, for example 



In the rest frame of the medium, we apply (|26 1 

S'^S' 



dEl, 



dt' 



gggCc At 
a 



IN) ' (H-4 ^ ^ I " > [>: 

R-S'^RS'^j C' x x 



(28) 

where the primed matrices, S', contain reflection and 
transmission coefficients evaluated at the rest frame fre- 
quency and wave- vector. 

The field amplitudes with different polarizations trans- 
form between laterally moving frames as follows (See ap- 
pendix A 2), 



a. 



M ■ a (IN) , 



(29) 



where, 



M 



1 



VyV 



V 2 V 2 V y —<T z 



(*1 

1 2 



has off-diagonal terms directly proportional to the mixing 
parameter 



V 



fen — OjVyky/c 2 



(30) 



The matrix, M is unitary, and the magnitude of the 
field intensity, a^a, is therefore only modified by the 
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FIG. 3: Contribution to the normal force (|36| from the mix- 
ing of polarizations due to the motion of the medium. Here 
we plot the quantity, (dP£ M) /dt) = £ £„ 2(-l)«» 7 (|72 OT | a - 



ir oc 



This represents the contribution to the normal force 



due to polarization mixing when a^f = ot^' , normalized in 
units of the incident momentum flux times V y /c. In figure (a) 
we plot {dP^ M) /dt) as a function of 6 for various 4>, where 
= (u/c)[cos(8)x + sm(0)(cos(4>)y + sin(0)*)]. In the 
main plot the radiation is incident onto a slab of thickness 
IOc/oj with e = 6.0 + O.Oli & fi = 1.0. The inset shows the 
case when the thickness is AQc/uj. Figure (b) illustrates the 
increase in magnitude of {dP^ /dt) when the refractive in- 
dex drops below unity: e = 0.1 + O.Oli & /j, = 1.0, where the 
medium becomes sensitive to the difference between s and 
p polarization at small angles of incidence (in this regime 
r\ — cot((9) sin(0) is large). Within both plots we used the 
reflection and transmission coefficients for a slab of dielectric 
of arbitrary thickness as can be found in |38| 



Doppler shift through the factor, (oj'/uj) 2 . Using (29 1, 



and Lorentz transforming ( 28 1 to the lab frame, we ar 



rive at the four-force experienced by a material in lateral 



motion, 



electric medium where S' 12 = S' 21 = is, 



dE M 
~dF 

dP M 

dt 




- S S^j tt(i N ) 

R-S^RS^ CxX 

U-s*s) Ciifeii 



where, A = A'/j, and, 

S = M^S'M 



(31) 



(32) 



The significance of (31) is that the four-force experi- 



enced by a laterally moving medium can be completely 
described by a unitary transformation of the rest frame 
scattering matrix S . The physics of radiation pressure 
on laterally moving media is entirely wrapped up in the 
dependence of the elements of S' on the Lorentz trans- 
formed frequency and wave-vector, generally referred 
here as to Doppler shift, and in the magneto electric ef- 
fect induced by the unitary transformation, M. Notice 
that when the angle of incidence is such that either k x or 
k z equal zero, then rj = 0, so that S = S\ and the only 
effect distinguishing the laterally moving medium from a 
stationary one is related to the Doppler shift [40] . 

Using the notation of (21 ), the unitary transformation 



/ 2 V 2 \ 

/ S^ + ^-^S^, v yV (s' 11 tr z -tT z S' 22 ) \ 



s = 



Vy-n ( a 'zS' 11 -S' 22 cr-) 



1+- 

„2,,2 



(33) 



/ 



There are several interesting features of ( 33 ) . Firstly the 
diagonal sub-matrices, Sn & S22 are only modified by 



the transformation ( 32 ) via terms that are second order 



in the velocity, and such terms are typically negligible. 
Secondly, the off diagonal terms linear in the velocity, 
S12 & S21, are related by S21 — cr z S\2(T z . Therefore, 
the mixed polarization reflection coefficients are related 
by a minus sign (72-21 = — ^12), while the corresponding 
transmission coefficents are equal {J~2\ — Tvi)- Assuming 
a rest frame scattering matrix which is left-right sym- 
metric, an examination of the off diagonal matrices in 
(33) shows that 7vL 2 i = —^21, ^12 = -^12, %.x = —T21, 
and 7i2 = — 7i2i which means that left-right symmetry 
still holds in this case as expected (recall the minus sign 
within the definition of the polarization shown in figure 

0- 



With the help of ( 33 1 we may now write out the four- 



force (31) for the simplest case of a plane wave incident 



of the scattering matrix given in (32 ) for a non magneto- on the slab from the left, 



dE M 
~dT 

dPu 

dt 



Ae c( x 



2 1 



Ae Q ( x 



2 1 



c 2 ) 



It' I 
1 1 ii 1 



+ (-1) 



gVyV 





\n> I 2 - 

1 /v 99 1 


It l 2X 

\'qq\ , 


Cxx\ 




\K' I 2 - 

1 /V 99 1 


It 1 

\'m\ 1 





+ (-1) 



qYvP- 



-a 



<+) 



(34) 



Comparing (34) with (27), it is clear that the lateral mo- 



tion affects the four-force (as discussed for the scattering 
matrix) through mixing the field amplitudes with differ- 
ent polarizations, and through Doppler shifting the fre- 
quency and wave-vector within the reflection and trans- 
mission coefficients themselves. 



It is instructive to limit our considerations to situa- 
tions of practical interest, i.e. for relatively small veloc- 
ities. Upon expanding (34) to first order in V y /c, the 



reflection and transmission coefficients are expanded as 



7?' ^7? 

qq qq 

work done, and normal force 



VykydlZ qq /duj. One then has for the rate of 



dE u 
~dt 



Ae c(x 



i-\n 



qq\ 



-\Tq 



>l'l\ 



2V y kyRe 



qq doj 



qq du 



I 2 



(-1)^ (|7^| 2 + \T qq ?) Re ( fl W^») (35) 



9 



'dP, 



dt 



Ae £ 



E 

{9} 



[1 + \K 



</'/ 



<IQ\ 



G 



(+)|2 
qL I 



T 



(-1)^(1 



\% 



I 2 ) Re (a;,, n 



qL a qL ) 



(36) 



In this limit we have, 77 = |fca;|fc z /A: 2 . The lateral force, 
proportional to the loss term in the round bracket on the 
right hand side of ( 34 1 , may be cast in the simple form 



dPu 

dt 



fell 
to 



dEu 
dt 



(37) 



Notice that within both ( 35 ) and ( 36 ) there are two veloc- 



ity dependent terms. One represents the change in the re- 
flection and transmission coefficients due to the Doppler 
shifted frequency response, while the other and less fa- 
miliar one is an interference term arising from the fact 
that the four force is sensitive to the phase difference be- 
tween the two complex polarization amplitudes a 1 ^ and 
a (+) 

This polarization dependent contribution is illustrated 
in figure[3]for a linear (electric) polarization as a function 
of the angle of incidence. The polarization dependent 
contribution to the force vanishes to first order in V y /c 
for circular polarization. In figure [3Jd we illustrate an in- 
stance where the contribution to the optical force due to 
polarization mixing may be enhanced. In this case we re- 
duce the permittivity to a value less than 1 and thereby 
make the medium sensitive to the difference between s 
and p polarizations for angles where rj = cot(#) sin(</>) is 
large. However, this method is certainly not unique, and 
it is likely that a greater enhancement would be possi- 
ble with something like the polarization sensitive mirrors 
of [26]. 

This is a peculiar effect indeed and does not appear 
to have been noticed before. When the phase difference 
is equal to an integer multiplied by it this contribution 
is maximized (either positive or negative), while a 7r/2 
phase difference causes this velocity dependent contri- 
bution to vanish. Furthermore, under special circum- 
stances the interference term and that due to the Doppler 
shift may cancel out so as to yield a velocity independent 
force. For a lossless medium this could occur e.g. when 
a[ + L } = a£\ and ck y Re[TZ qq dTZ qq / 'du] = (~l) q r,\TZ qq \ 2 . 

The above findings have a relevance to the recent 
controversy regarding the existence of a drag component 
to the Casimir force, discussed in the introduction. 
Firstly, if the field has an ill defined phase (e.g. a 
thermal field, or the quantum vacuum), then the first 
order effect of polarization mixing should average to 
zero, while the Doppler shift will not, c.f. [S]. Secondly, 
the contribution of polarization mixing to the force 



vanishes when the medium has a degenerate response to 
the two polarizations, whereas the contribution due to 
the Doppler shift does not. Finally, it is worth pointing 
out that the lateral force due to first order polarization 
mixing, and that due to the Doppler shift have a differ- 
ent dependence upon the direction of incidence, and so 
would not be expected to cancel in a general situation. 
From ( 35 1 and |37|) it can be seen that the velocity 



dependent force in the y direction has a contribution odd 
in k y coming from the term proportional to rj, and even 
in k y from the term proportional to the derivatives of 
the reflection and transmission coefficients. An isotropic 
field will therefore — through the effect of the Doppler 
shift alone — exert a lateral force on a body that is in 
motion, so long as that body can absorb photons. This 
is related to the phenomenon of universal drag that 
was highlighted a few years ago [27], where a dielectric 
in a thermal field was generally found to experience a 
frictional force that is linear in the imaginary part of the 
dielectric susceptibility. 



B. Motion along the x axis 

When a plane-wave is incident onto a surface that 
moves with velocity V — V x x parallel to the surface nor- 
mal, we cannot in general use a monochromatic transfer 
or scattering matrix to determine relationships between 
field amplitudes. A plane-wave of given frequency in the 
lab frame is reflected with a different frequency. In the 
slab rest frame the amplitudes of waves moving to the 
left (o^r') and moving to the right (a^t') suffer distinct 
Doppler shifts. Therefore for plane waves of fixed fre- 
quency u> in the lab frame impinging on both sides of the 
slab, there are eight rather than four output amplitudes 
at three different frequencies (w,oj (+, ,u; <-) ). The result 
( 26 1 no longer applies and must be modified (see figure 



To apply (26) to a slab in motion along x, we note 



that, due to the time averaging of sectionpl] that leads to 
([3]), the two different input frequencies observed in rest 
frame of the slab contribute additively to the energy- 
momentum tensor, and hence the force |29j . Therefore, 
we sum the force due to incidence from the right and left 
of the slab separately, 
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dt' 



dtf 



A'e c , f 



r<±) 



{±> 

^ / tt^' f J2 - Spies') 7r<±' (' x ^x 

Vc' (±) v / 

^ \ tt^ (l 4 - S'+S' ) tt (±) Cn (±) fe|| 



(38) 



The components of the rest frame scattering matrix, S , 
are evaluated at the two different rest frame frequencies, 
u/ (±) = 7 (w =F V^lfc^l) depending on whether they are as- 
sociated with incidence from the left (+) or from the right 
(-). For the reflection coefficients of r -polarized wave im- 
pinging from the left one has e.g. 1Z' qr — lZ qr (uj' {+y ) and 



1Z' = lZ qr (ijj'( } ) for waves impinging from the right. 



A Lorentz transformation of ( |38| into the laboratory 
frame gives the dependence of the observed four-force on 
the velocity of the medium. The rate of change of energy 
and normal force on the medium are given by (A! = A), 



dE M \ _ ^oc a n £(±M±) 



dt 



l IN / j Sj; 

{±} 



1 4 + C 



(39) 



dt 



Ae 
2 



/t ^£(±)„.(±) 
{±} 



«C5 



'(±) 



r,. 



1j s' 



(40) 



and the lateral force is given by the same expression as in 
the rest frame ( 38 ) . All that remains is to write all quan- 



tities in terms of those in the lab frame. The directional 
terms £, in particular, transform as, 



'(±) 



c 



c|A4 (±) 



'(±) 




w /c±) 7 (w=fVx|*x|) 



In this case the field amplitudes of different polarization 
do not mix between reference frames, and are related by 
the diagonal matrix (see appendix All, 



ad = 7 Ua - 



a„ 



As in the previous section we examine the simpler case 
of incidence from the le ft onto a medium where S' 12 = 
S' 21 = 0, so that |39||40| become 



dE u 
dt 

dP M 

dt 



Acqc ( coj'k' 



Aeo ( coj'k'z 
~2~ 



{<?} 



{<?} 



1 - \K' 



99 I 



^0 



It' I 

I '99 I 



qq 1 



7W 



'0 



1 + 1^:. 



</</ 1 



IT I 

I 'qq I 



2 



1^: 



<i<i 



2 _ | T , |2 

I 'qq I 



99 I 



T 
'99 1 



(41) 



where the superscript (+) on the rest frame quantities 
has been omitted. The velocity dependence of the force 
in (41 ) has contributions coming from the rest frame loss 



term, 1 — \TZ' \ — \T qq \ , the dispersion of the reflection 
coefficients, and a prefactor of cuj'k x /uj 2 that is < 1 for 
V x > 0, and > 1 for V x < 0. It is important to note 
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that the force ( 41 ) comprises a lateral component that 

I 1 2 i 1 2 

depends on the loss term 1 — \TV \ — \T qq \ ■ For a slab 
moving along x, however, such a loss term is clearly not 
the only contribution to the rate of loss of energy from 
the field and — at variance with the case of a laterally 



moving slab — the lateral force component in (41) cannot 



be cast in the simple form (37) 



Retaining only terms linear in V x /c in (41 1 and using 
the notation of section III A[ the rate of change of energy 
and normal force on the slab is, 



lit 



Ae c ■ 



{«} 



c x (i~\n qq \ 2 ~\T qq \ 2 ) 



2( x V x k x Re[K qq ^ 
V x 



dT* 
T u 'qq 

lqq duj 



l-(l + 2( 2 x )\K qq \ 2 -\T q 



'991 



„(+)|2 



(42) 



dpil 

~dt 



Ae Q\k\\ 



E 

{9> 



{1 - \K qq \ 2 - \T qq \ 2 ) Cx- 



dU* dT* 
2CVkHelTZ qq + T qq 

*>*X '1%™ I l^-qq ' lqq 



„(+)|2 



(43) 



dP ± 



Ae ( x x 



E 

{9} 



dlZ* dT* 

Cx (1 + \K qq \ 2 - \T qq \ 2 ) - 2Q x V x k x Re ( K aq ^- - % ! " 



^ duj ,qq Olo 

v x 



c 



l + 3\K q q\-\T q , 



„(+)|2 



(44) 



There are several comments to be made about (42 
44 ) . Firstly, there are no interference terms between TE 



and TM polarizations due to the motion, as there were in 
( 35 -[36 ) , and so the phase of the incident field is irrelevant 
to the optical force. Secondly, we have the term propor- 
tional to — (l + 3\lZ qq \ 2 — \T qq \ 2 ) in the normal component 
of the force. This quantity is always negative, and so the 
term acts like a friction against the acceleration of the 
material by the field. Indeed, for normal incidence onto 
a dispersionless, lossless medium, this contribution be- 
comes — AV x \R, qq \ 2 / c, which is the radiation friction (or 
damping) term investigated within |17l 118] . The contri- 
bution of the dispersion of the reflection and transmission 
coefficients within (42 44) is not fundamentally differ- 
ent from that found for lateral motion, and allows for 
a damping or amplification of the centre of mass motion 
depending upon sign of the gradients of the reflection and 
transmission coefficients [20l [39] . Note that, unlike the 
case of lateral motion where the polarization mixing pa- 
rameter r\ depends on the plane of incidence, here the role 
of the plane of incidence of the field is not fundamentally 
altered from a medium at rest. 



IV. CONCLUSIONS 

We have developed a general theory of radiation pres- 
sure for planar media in terms of a 4 x 4 scattering ma- 
trix, which allows for coupling between the s and p po- 




CO 



FIG. 4: A slab in motion along the x axis does not mix polar- 
izations, but mixes frequencies. Incidence from the left, when 



(-) 



<UJ, 



V x > yields a reflected amplitude at a frequency, ui 
while incidence from the right yields a reflected amplitude at 



,(+) 



> lu, which breaks reciprocity. 



larizations in transmission through and reflection from a 
medium. Applying this theory to the problem of the ra- 
diation pressure experienced by a planar medium in mo- 
tion, it was found that for lateral motion (see figure [2]) 
the effective scattering matrix of the medium is such that 
polarizations are coupled (sec (33)), and something akin 



to the 4x4 scattering matrix description is unavoidable. 
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This has unusual consequences for the optical four-force, 
which then has a velocity dependence that changes ac- 
cording to the phase difference between the amplitudes 
of the two different polarizations (see (35 37)). Indeed 
in the general case (31) the dependence of the force on 



the rest frame field amplitudes is quite involved. Mean- 
while, for motion parallel to the surface normal there is 
no mixing of the polarizations, and to the lowest order in 
V/c our formalism recovers the known results for moving 
surfaces [17l HH [20] . It is indeed possible to apply this 
formalism to the general case of a medium in motion in 
an arbitrary direction, for we can consider this as a com- 
bination of motion normal and parallel to the plane of 
incidence. 

One clear advantage of this formalism is its relative 
simplicity. Work on the problem of quantum friction 
often encounters very cumbersome formulae on which 
there seems to be no general agreement (e.g. [H]). The 
fact that lateral motion can be thought of as hav- 
ing two effects — a unitary transformation of the trans- 
fer/scattering matrix, and a Doppler shift of the fre- 
quency within the reflection & transmission coefficients 
is therefore advantageous, and it should be possible to 
apply this formalism to the problem in a straightforward 
way through a generalization of [31 . 
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we obtain the equations 



s '(±) '(±) 



-2.r 



a 9 



2x 



(In 



(A2) 



and 



~'(±) '(±) 



;(±) 



s '(±)„/(±) 

2 



+ e 2y 'a 



V x 



~2z 



(±) \ „,(±) 



;(±) 

'2y 



„(±) 



(A3) 



with the third of ( Al ) holding identically, given ( A2 - A3 



We then insert the components of the polarization unit 

(±) - 



vectors which from (|6||7| are; e^jr = ck /uj; e{ 



k y/ k \\' 4? = Tc\k x \k z /u}k\\; eff 



-k z /k\\ ; and e^ 1 



^c\k x \k y / Lok\\ . For (A2) this gives 

,(±) 



'(±) 



(±) 



(A4) 



where L)^' — ^{lu =f | /c^ | ) . Note that we do not con- 
sider the case where k x changes sign between reference 
frames, although this is interesting. Inserting (A4) into 



( A3 ) along with the components of the unit polarization 



vectors, we find 



,/(±) 



,(±) 



„(±) 



The field amplitudes consequently transform as 

V x \k x \ 



7 1 



(A5) 



(A6) 



Appendix A: The Lorentz transformation of 
polarization 

Here we examine how polarization transforms between 
reference frames. We consider motion of two kinds; ei- 
ther parallel to the surface normal (motion along £), or 
orthogonal to the surface normal (motion along y). 



between transversely moving inertial frames. In this case 
the polarization of the wave does not change between 
frames, the amplitude being either increased or decreased 
according to the ratio of the rest frame and the laboratory 
frame frequency. 



2. Motion along y 



1. Motion along x 

When a planar medium moves in a direction paral- 
lel to the surface normal, then after scattering from the 
medium, the fields Q will no longer be monochromatic. 
However, the phase factors within Q will take the same 
value in any inertial reference frame. Therefore after ap- 
plying the transformation 



E ' x - E x 

E'y = l(Ey - V X B Z ) 

E' z = j(E z + V x B y ) 



(Al) 



In the case of a medium moving laterally as in figure 
[2j the frequency to is conserved. Applying the transfor- 
mation 



E' X = 1 {E X + V V B Z ) 



E 'y - E v 



E' Z = 7 (E Z -V V B X ) 



(A7) 



to Q along with the equivalent primed quantities leads 
to the equations 



s /(±) /(±) _ 



(±) / £ (±) 



V, 



» £ (±) 



V, 



(A8) 
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and 



~'(±) '(±) 



*/(±) /(±) -(±) (±) 

°2y a 2 = e ly <*1 



The third equation in (A7) holds identically when we 



apply (A8-A9). We insert the components of the unit 
polarization vectors given in appendix |A 1| In the case 
of ( A8 1 this gives 



a 



'(±) 



70/ 



„(±) 
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^)T^\k x \k z a{ ±} 



wAi||fcj 

Defining 77 = \k x \k z /(k 2 — VykyUj/c 2 ) and noting that 

fell fen 



- Vykyljj/C 2 ) 



we find the transformation formula for the p-polarized 
amplitudes 



'(±) 




(±) _,_ v yn (±) 



^l + Vyhf/C 2 



(A10) 



Inserting (A10) and the components of the unit vectors 



into ( A9 1 then gives the transformation formula for s- 



(A9) polarized amplitudes 



'(±) 




yjl + V 2 V 2 /c 2 



In matrix notation, (A10) and (All I become 



-IM.q 

Ul 



where 

M = 



1 2 



VyV 



a, 

la 



(All) 



(A12) 



(A13) 



This unitary matrix, M describes the transformation of 
the polarization of a plane wave between laterally moving 
reference frames. 
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